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We study the spectrum, the massless S'-matrices and the ground-state energy of the 
flows between successive minimal models of conformal field theory, and within the sine- 
Gordon model with imaginary coefficient of the cosine term (related to the minimal models 
by "truncation"). 

For the minimal models, we find exact ^-matrices which describe the scattering of 
massless kinks, and show using the thermodynamic Bethe ansatz that the resulting non- 
perturbative c-function (defined by the Casimir energy on a cylinder) flows appropriately 
between the two theories, as conjectured earlier. 

For the non-unitary sine-Gordon model, we find unusual behavior. For the range of 
couplings we can study analytically, the natural ^-matrix deduced from the minimal one 
by "undoing" the quantum-group truncation does not reproduce the proper c-function 
with the TBA. It does, however, describe the correct properties of the model in a magnetic 
field. 



April 1993 



1. Introduction 

One of the most useful methods of understanding an integrable two-dimensional field 
theory off the critical point has been to find the exact S-matrix. Strict consistency cri- 
teria have allowed exact S-matrices for a large number of massive integrable models to 
be conjectured, using intuition from perturbed conformal field theory. Many have been 
essentially verified by a number of checks, especially that of the thermodynamic Bethe 
ansatz, where one can calculate the Casimir energy from the S'-matrix and check that it 
is correct at the conformal point. 

Massive theories are not the only type of field theory which can be studied by these 
methods. There are numerous situations where perturbation of a conformal theory by a 
relevant operator causes it to flow to another conformal theory. In the midst of the flow, 
the model has massless excitations (the only ones to survive in the infrared limit) but is 
not scale-invariant. Thus one can hope to find an exact S'-matrix for these excitations, and 
calculate the associated Casimir energy in order to study the flow. This has been done for 
several models in |l]||. There are subtleties in defining an S-matrix for massless particles, 
but since our models are integrable and have an infinite number of conserved currents, we 
expect that the idea of massless-particle scattering can be implemented. In particular, it 
was shown some time ago that in the lattice XXX model (which is at a conformal point), 
one can define an S-matrix by deriving phase shifts of quasi-particle excitations from the 
Bethe ansatz ||. 

One crucial point to realize about a massless S-matrix is that it should be viewed 
as describing an (irrelevant) perturbation of the infrared fixed point. The reason for 
this is simple: the massless excitations are those which remain in the spectrum at the 
infrared fixed point. All known massless S-matrices follow this pattern; one mainly uses 
intuition gained by studying the infrared, not the ultraviolet, conformal field theory. We 
can even make this observation precise. In these massless flows, the mass M provides the 
scale; M = gives the ultraviolet fixed point while M — ► oo gives the infared one. Since 
the particles are massless, they are either right- or left-moving; we parametrize the left- 
movers' momenta by pl = — (M/2) exp(— 9), and the right-movers by pr = (M/2)exp#. A 
Lorentz-invariant S-matrix element Sll describing scattering of two left movers depends 
only on the ratio of the two momenta, so it depends only on 9± — 62 and not on M. 
We can always rescale M — > 00 by shifting the rapidities. The LL and RR S-matrices 
are independent of this shift (although Slr is not), so they are characterized solely by 



1 



properties of the infrared fixed point. In fact, one can think of the LL and RR S-matrices 
as being the S-matrices for the conformal field theory. The LR S-matrix does depend on 
M, so this is affected by the perturbation. 

There are two purposes of this paper. The first is to describe the flows between 
the minimal models of conformal field theories, completing a picture initiated in 
We will show that the spectrum consists of massless kinks, which can be thought of as 
interpolating between the vacua of a scalar field with a degenerate- well potential. Their 
S'-matrix is related to a well-known RSOS one. We will verify that these S-matrices 
give the appropriate ground-state energy (on a cylinder) in the IR limit in all cases, and 
for all scales in the simplest cases. The ground-state energy provides a non-perturbative 
c-function which allows one to see explicitly the crossover from one minimal model to 
another, a phenomena which had been earlier studied perturbatively 

The second purpose is to discuss the related S'-matrix simply obtained by "undoing" 
the quantum-group truncation, which may describe scattering in the sine-Gordon model 
with imaginary coefficient of the cosine term. Besides its relation with minimal models, 
this theory is interesting in its own right: it is the continuum limit of the 0(n) model in 
the low temperature phase, and for n — > 0, it describes the "dense" phase of self-avoiding 
polymers. As we have discussed at length in our previous paper 0, this non-unitary 
sine-Gordon model has a number of interesting but difficult-to-handle properties. Our 
discussion will further elucidate these properties, although we will see that the S-matrix 
does not seem to completely describe the model. It does reproduce the correct behaviour 
with a magnetic field in infinite volume, but we will show (at least in the simplest cases) 
that for properties on a cylinder (i.e. at finite temperature) the behavior in the ultraviolet 
is incorrect. 

Various scenarios to explain this feature are proposed. 

2. The flows between the minimal models 

2.1. S-matrix 

The models we study in this section are the minimal models of conformal field theory 
perturbed by their least-relevant operator $13: 

S = S* min (t) + 8/3 J $ 13 d 2 x, (2.1) 
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where the unperturbed conformal field theory has central charge c = 1 — 6/t(t + 1). The 
physical properties depend crucially on the sign of 6(3. For 6(3 > 0, the model is massive. 
The underlying lattice model || has t degenerate vacua |§. An effective Landau- Ginzburg 
theory |K| describes the theory by a single scalar field with a potential with t degenerate 
wells. The particles are kinks which interpolate between adjacent vacua, and an S-matrix 
consistent with the quantum-group symmetry was conjectured in |Tl|,12|. This S-matrix 
was checked by showing that the resulting exact Casimir energy gives the correct result at 
the conformal point PJl^l- 

For 6(3 < the model is qualitatively different: it flows from the t-th minimal confor- 
mal field theory to the (t — l)-st one HH. In the midst of the flow, the model should have 
massless excitations. In the t = 4 case (the flow from the tricritical Ising model to the Ising 
model) the spectrum was shown |I| to consist of a left-moving and a right-moving particle, 
which can be thought of as the massless Goldstinos arising from spontaneously-broken 
supersymmetry. Their S-matrix was found, and the ground-state energy was calculated. 
The effective central charge was shown to smoothly flow from 7/10 in the ultraviolet to 
1/2 in the infrared. In addition, the ground-state energy was also conjectured for all t 
without benefit of an S'-matrix. This was not omniscience; looking at the massive case for 
all t and the massless case for t = 4 makes the guess straightforward. This was shown in 
to have the correct properties. 



Obviously, it would be useful to find the S'-matrix for all t, and to then verify that 
the ground-state energy conjectured in JI] is correct. We will find the S-matrix for all t, 
showing that it agrees with the already-known t = 4 case. We find that in the general case, 
the spectrum consists of kinks reminiscent of the massive case, but here they are massless. 
Although one is used to thinking of kinks as field configurations with non-zero energy, there 
is no reason why they cannot be degenerate with the ground state in a quantum theory. 
We will explicitly calculate the ground-state energy for t = 5 and find the conjectured 
result. Technical complications prevent us from doing so for all t, but we show that in the 
infrared limit our S-matrices give the correct answer. 

We must first find the appropriate particle content. For guidance, we can again study 
the underlying lattice model ||. As discussed in M, whereas the massive phase (regime 
III) has t — 1 degenerate vacua, the critical point described by the t minimal conformal 
field theory separates this from a phase with t — 2 degenerate vacua (regime IV). The 
result can be seen in the effective Landau-Ginzburg picture: for even t, the central vacuum 
disappears, and for odd t the two central vacua coalesce into one. The dimension of the 



3 



perturbing operator is the same in both phases, and since there is only one such operator 
in the conformal theory (the $13 operator), the two phases must correspond to perturbing 
by $13 with different signs. So the clear guess is that as in the massive case, there are 
kinks interpolating between adjacent vacua, but for a given t there is one less vacuum. 
In addition, the kinks here are massless, so for every pair of adjacent vacua we have a 
right-moving and a left-moving kink. 

We derive the simplest S-matrix for this spectrum. Since scattering is diagonal in the 
L and R labels, we can consider the two-particle S-matrices Slr, Sr l and S l l = Srr 
separately.0 The mass scale M oc \5 (3\ <yt+1 ^ ^ . The S'-matrix depends only on the rapidity 
difference 9 = 9\ — 62- for LL scattering, the ratio of the two left-moving particles' mo- 
menta is the only Lorentz-invariant kinematic quantity, and for RL scattering the Lorentz 
invariant s = (p\ +P2) 2 also depends only on the rapidity difference. In order that the scat- 
tering be factorizable, the S'-matrix must obey Yang-Baxter equation for a three-particle 
state: 

$12(61 — #2) £13(01 — 63)823(62 — #3) = £23(^2 — 63)813(61 — 6s)Si 2 (6i — 62) (2.2) 



where S12 is the two-particle S- matrix acting on particles 1 and 2. A two-kink configuration 
can be labeled by three vacua; a two-particle S-matrix element can be labeled by four 
because only the middle vacuum can change in a collision. The simplest possibility for 



an LL or RR S-matrix obeying ( |2.2| ) and having this structure is the well-known RSOS 
S-matrix [Hi 
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sinh / u(zm7r — 9) 



(2.3) 



1 It may appear absurd to define an S'-matrix element Slr for a right-moving particle starting 
to the right of a left-moving one. However, it should be interpreted as a matching condition on 



plane-wave eigenstates of the Hamiltonian [15]. 
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where 

(3 m = sinh(zm//7r) 

and |U = l/(t — 1). The horizontal line denotes a right-moving particle with rapidity 
#i, while the vertical (moving from bottom to top) has rapidity 62, and 6 = 61 — 62- The 
allowed vacua run from 1 to t — 2, and notice that as required, adjacent vacua differ only by 
±1. The factors of (/3 m ) l6 '/ 7r ensure crossing symmetry, but can be absorbed by a rapidity- 
dependent change of basis. In the massive high-temperature phase the scattering is given 
by the RSOS .S-matrix (|2.3|) but with \x = 1/t |TI,12|. This of course is a manifestation 



of the fact that there are t — 2 vacua in the low-temperature phase, but t — 1 at high 
temperature. 

In p| , two types of RL S'-matrices were discussed: one proportional to the identity and 
the other Sll- The same possibilities can be considered in our case, since both trivially 
still obey (|2.2|) , but they are not correct here (the TBA gives the wrong central charge 
in the UV limit), probably because they have too much symmetry. There is a "restricted 



quantum affine symmetry" [[[6[ which does exist in the critical limit but which is broken 
in the midst of the flow. Consider for instance the t = 5 case and the flow from the three- 
state Potts model to the tricritical Ising model. The restricted quantum affine symmetry 
of the infrared fixed point, the tricritical Ising model, is in fact iV=l supersymmetry. The 
flow approaches this fixed point by the irrelevant operator $3,1, which in this case is the 
supersymmetry operator, of dimension (3/2,3/2). This operator is the lower component of 
a superfield with the energy-momentum tensor. In order to perturb a model and preserve 
the supersymmetry, one can only perturb by the upper component of a superfield. Hence 
in this case the supersymmetry is explicitly broken off the critical point. As explained in 
the introduction, the LL and RR S- matrices essentially describe only the infrared limit, so 
Sll and Srr should have this symmetry (as they do: see [^TJ,|T2|), but Srl should break 
it. This excludes immediately the diagonal choice for Srl- The second choice Srl oc Sll 
was proposed in JT7J where it was argued to preserve "diagonal supersymmetry". We are 



not sure what such a concept means, but we have checked explicitly using the TBA that 
the corresponding UV central charge is incorrect (we find cjjv = 9/5). Similar arguments 
apply for all values of t by replacing N = 1 supersymmetry by the appropriate restricted 
quantum affine symmetry. For the second choice, also proposed in [fT7| , we find indeed 



incorrect values of cuv (except for t = 4 of course). The two choices, Srl proportional to 
the identity or to Sll, therefore must be excluded.i 

There is fortunately a whole family of possible RL S'-matrices which obey the Yang- 
Baxter equation ( |2.2| ). Consider the situation when one of the three particles in (|2.2 



is 



right-moving while the other two are left-moving. Here one of the S'-matrix elements in 
( |2.2|) is of LL type, while the other two are of RL or LR type. If we shift the right-mover's 
rapidity from 6i to $i — ia, ( |2.2| ) is still of course obeyed. The key observation is that the 
S'-matrix elements affected by this shift are of RL type if +9i appears in the argument 
and of LR type if — 0$ appears. The equation is then equivalent to the all-LL problem, 
and hence is solved, if 

Srl oc S ll {9 + ia) S LR oc S LL (9 - ia). 

To fix a, we make the simple demand that Slr = Srl, because we know that in the 
t = 4 case this holds JIJ. Obviously, the already-ruled-out a = would work, but there is 
another interesting possibility. This is to set a = tt /2/x. We are also free to multiply by 
an overall function, so we set 



This replaces Z with Z and sinh with i cosh in each of (^.3|) . 



2.2. "Undoing" the quantum-group truncation 

We could find Z(9) and Z{6) directly by demanding the usual constraints of unitarity 
and crossing symmetry. However, we first find them for a related, simpler S'-matrix. A well- 
known "reduction" procedure then gives the RSOS .S-matrix from this S'-matrix [|T9| , |12"|| . In 
the context of integrable lattice models, this is known as the vertex-IRF correspondence. 
The unreduced S'-matrix we find will also be used in the next section, where we study the 
sine-Gordon model when the coefficient of the potential imaginary. For the moment one 
can think of it as an intermediate object that makes the study of normalization convenient. 



2 In fact, the Srl proposed in [fTf! for the flows into the SU(2)k (k > 1) WZW models do not 
give the correct cuv = 3. The proposed Sll and Srr are correct (as shown previously in [[l8f ). 
but Srl is equal to the one discussed below. 
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The unreduced ^-matrix describes the scattering of two doublets, one left-moving and 
the other right-moving. For the LL or RR scattering, the only possible ^-matrix for a 
doublet (u, d) consistent with factorizability (|2.2|) is the usual sine-Gordon S'-matrix [|0 



a = Z(9) sinh/x(6> — in) 

b = Z(d)smh/j,d. (2.5) 
c = —Z(9) sinhz//7r 

where a describes the process uu — > uu, b describes ud — > ud, c describes the reflection 

ud — > du and there is a symmetry under u <-> d. Again, [i = l/(t — 1). Unitarity and 
crossing fix Z to be 

1 i r°° dk sinh - 1) 

Z{9) = - , n rexp- / ^ S mk9 . 2 V , (2.6) 

v ' smh(j,(9-nr) 2j_ QO k sinh g cosh ^ V ; 

This generalizes to the XXZ case the massless LL ^-matrix for the XXX model (/x — > 0) 
. It also arose in the WZW models of @] . 

The unreduced Srl then follows from ( |2.4j ). (This S'-matrix has the slightly odd 
property of having crl = —clr-) To fix Z(6), we use the unitarity and crossing relations 
0,| . These are 

a LR (0)a RL (6) = b LR (9)b RL (9) + c LR (9)c RL (9) = 1 
b L R(9)c RL (9) + c LR (0)b RL (9) = 

(2.7) 

b RL {9)a RL {m + 9)+ c RL {9)c RL {m + 9) = 1 
b RL {9)c RL (m + 9) + c RL (9)a RL (in + 9) = 0, 

where the first two apply only for 9 real, and the last two also hold true for LR-type 
elements. The second and fourth are satisfied for any a and Z. The other two give 

Z(9)Z*(9) = — -— * — = Z(9)Z(in + 9) (2.8) 

cosh ^(ztt - 9) cosh /j,(9 + in) w v ; v ; 

The simplest solution of these is 

^ ^ cosh/x(6> — in) 
fr T 2 (^ + (2n - 1)^ + g)r(l + (2n - 2)^ - g)T(j + 2n^ - g) 
Al r2( i + (2 „ - 1)/. - g)r(± + (2n - % + g)r(± + 2n^ + g) ( 2 ' 9 ) 

1 dfc sinh^ 



To find the appropriate RSOS S-matrix, we rewrite the LL S-matrix as 

S LL = Z(9) [sinh^(# - m)P + sinh(//0)e] , (2.10) 
where P is the permutation matrix, and e is 
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The S-matrix ( |2.10| ) is equivalent to ( j2.5| ), but we have redefined the states \u(8i)) — > 
e~ e ^ 2 \u(9i)) and \d(9i)) — > e e ^ 2 \d(9i)). This ensures that e obeys the Temperley-Lieb 
algebra |H| . Any S- matrix of the form ( |2.10|) satisfies ( |2.2| ) as long as e obeys this algebra. 
Our RSOS S'-matrix is given by fl2.10|) , but with an e which satisfies the same algebra but 
acts on the RSOS particles. The elements of this new e are 

[sin(a/i7r) sin(c / U7r)] 1 ^ 2 



b 


c 
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d 



-5 



bd- 



sin(6/x7r) 



where adjacent vacua must differ by ±1. This gives (|27|) after we do another rapidity- 
dependent "gauge transformation" to get the (/3 m y 6 / n and the minus sign. Slr = Srl 
follows from (|2~4[) . 

We have completely determined the RSOS S'-matrix, because we now know Z(9) from 
( |2.6| ) and Z(0) from (|2.9| ). We can check that it agrees with the previously-known case 
t = 4 (/i = 1/3) |l|. In this case, the RSOS structure is trivial: all one can do is go back 
and forth between the two vacua. Thus the spectrum consists of one left-moving and one 
right-moving particle: the Goldstinos of 0. The two S- matrix elements Sll and Srl 
follow from either of the last two relations in ( |2.3| ). To show that this S- matrix agrees 
with that of 0, one uses the relation 

sinMag+jTr) . /"» dk z . nk0 ^m- 2 \ (2n) 

smb. i_i(a9 - in) 1 J_ OQ k sinh ^ L 

which is good for a > and < fi < h. This yields Sll = 1 and Srl = — tanh(6*/2 — iir/A), 
as required. 

To check that these S-matrices indeed describe the scattering for all the flows between 
minimal models, we will calculate the ground-state energy E(R) on a circle of radius R. At 
the ultraviolet and infrared fixed points this is proportional to the central charge c of the 
appropriate field theory by using the relation E(R) = ttc/6R. This calculation was done 
in for the case t = 4; in addition, the appropriate equations for all t were conjectured 
there. We will derive the equations of M explicitly for t = 5. 
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2.3. TBA 



The technique we will use to calculate the ground-state energies from the S'-matrix is 
known as the thermodynamic Bethe ansatz (TBA) ||15|| . The TBA gives the free energy 
/ = —InZ/R of a particle gas on a line of large length T at a temperature 1/R. This is 
equivalent to considering a spacetime of a cylinder of radius R and length T. If we change 
our point-of-view and think of the spatial direction as (Euclidean) time, only the ground 
state contributes to Z, since only the lowest-energy state can propagate over the large 
distance T. Thus Z = exp(—TE(R)). The equivalence between the two approaches (often 
called modular invariance) means that E(R) = — In Z/T = Rf /T. 

To calculate the free energy, we need to diagonalize the "transfer-matrix", meaning 
that we find the eigenvalues for scattering one particle through an ensemble of all the 
others. Since the t = 4 massive case and our t = 5 massless case are very similar, we first 
recap the results for the t = 4 massive case ||]. There are three degenerate wells labeled 
(+, and — ) with a massive kink interpolating between each adjacent pair. The S'-matrix 
elements are given by (|2.3| ) with \i = 1/4. The transfer matrix T(9\9i, 02, ■ ■ ■ 9n) can be 
represented pictorially by 






± 





± 





± 





± 





± 



6>i O2 On 

where the intersections represent S'-matrix elements, and the ± denote vacua which can 
be either + or — . We impose periodic boundary conditions, so N must be even and the 
rightmost + or — is the same as the leftmost. This is equvalent to the lattice Ising model 
transfer matrix: the configurations are specified by the ±, and the S'-matrix elements are 
the same as Ising-model Boltzmann weights. 

To find the transfer matrix eigenvalues A(6>|6>i, O2 ■ ■ ■ On), we follow |2j] and Q. We 
can ignore the (f3 m ) ±ze ^ n pieces, because they cancel out in T(6\0i, #2, ... , On)- Using the 
weights ( p.3|) one finds that if one ignores the Z{9), shifting 9^9 + 2iit is equivalent to 
flipping (+ «-> — ) one of the two Ising spins in a weight and multiplying the weight by i. 
In the transfer matrix this corresponds to flipping, say, the bottom row of spins. Thus 

H0\m) - ^P^^=r(-l)^X(9 + ^\{9 t }l (2.12) 
1 lz=i Z\9 - 9i) 
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where r is the eigenvalue (±1) of the configuration under spin-flip. We note that 
exp(iV6> / '4) \(9\{9i}) is a polynomial of order N in exp(6>/2), and by the relation fl2.12p , the 
odd terms vanish for r = (— 1) N / 2 and the even ones for r = -(-l) N / 2 . Thus A(0|{0i}) 
must have N/2 zeros in the strip —n < Im9 < n in the first case and N/2 — 1 in the 
second. By analyticity, it can thus be written as 



where A is independent of 9. Specifying the zeros {xi} (which depend on the rapidities 
9\, 92 ■ ■ ■ 9n) specifies the eigenvalue. To find the zeros, we use the "inversion" relation i 
derived in plfl : 



9-9,, " . . ,9-9, 



A(0|R})A(0 + z 7 r|{0 l }) = ncosh(^)+rnsmh(^-^). (2.14) 



i=l i=l 

r N 



Thus if x is a zero it must satisfy Yli=i tanh((x — 9i)/2) = —r{—l) N l 2 . The N solutions 
to this polynomial equation in the strip are of the form Xi = yi ± ift/2, where yi is real. 



The N/2 (or N/2 — 1) zeros {x,} needed to find an eigenvalue X(9\{9{}) in (|2.13| ) are then 
specified by choosing one from each pair. (If we choose both the left-hand-side of ( |2.14| ) 
will not be zero for some other x,, whereas the right-hand side still will.) 

We have found the eigenvalues of the transfer matrix for the t = 4 massive case, and 
one can then calculate the ground-state energy as in || . Generalizing this to our massless 
t = 5 case is only slightly more complicated. The transfer matrix displayed in the diagram 
above must carry an additional label corresponding to L or R on each line. The parameter 



H = 1/4 in both cases so Sll = Srr is the same, and Srl is given by the shifted fl2.4T) . 
Thus, just as we did in the Yang-Baxter case, we can reduce this to the LL-only case (the 
problem just solved) by pulling out the Z{9) as in ( |2.12| ) and shifting the rapidities of all the 
right-movers by 2m. This shift affects almost nothing of the above diagonalization. The 
eigenvalue \(9\{9i}) must have the same number of zeros in the strip, so ( |2.13| ) is exactly 



3 This relation is not hard to derive: one first finds T(0)T(6 + itt) explicitly for N = 2 but 
without periodic boundary conditions. One then builds up the relation for any N with periodic 
boundary conditions. 
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the same, and the only effect on ( |2.14j ) is to multiply the right-hand-side by (—1)^. The 
eigenvalue zeros Xi = yi ± in/2 are therefore the solutions of 

ji sinh(!^ - % ) H mM*^ - f ) 

With this diagonalization complete, we can find relations for the level densities and 
particle densities. The total numbers of left and right particles are each conserved, so we 
define the densities of states Pl(9) and Pr{9). In addition, we define the density P(y) 
of solutions to ( |2.15| ). A given particle configuration (eigenstate of the transfer matrix) is 
specified by the distributions of rapidites for L and i?, and the distributions of eigenvalue 
zeros for Xi = yi + in/2 and for Xi = yi — in/2: we thus define the associated densities 
Pl(0), Pr(9), P+(y) and P-(y). As we have shown, P(y) = P-(y) + P+(y). Taking the 
derivative of the logarithm of ( |2.15| ) gives 

2nP{y) = [ d9 PlW + P«V\ (2.16) 
yy> J cosh(y-6) v ; 

To find Pl{9) and Pr(9), we quantize the momenta: the periodic boundary conditions 
require that for an eigenstate of the transfer matrix 

e i * il A(6 i \6 u 62,...,6 N ) = l. (2.17) 



Taking the derivative of the logarithm of (|2.17 ) and using (|2.13 ) gives 



2nP L {9) = ^le- e +-^-Im In A(0|0i, ...,9 N ) 
2 au 



jle-°+ J M'W - e')p L (0') + f M'W - 0')pr(0') (2.18) 



where 



cosh(y — 9) 

(f)(9) =Im^-\ivZ{9) = - [ dk cos k9 

^ ' d9 1 ' 8 J cosh 2 f 

~ m =Im^\nZ{9) = m, 



(2.19) 



where we have used the relation ( |2.11|) . The explicit expressions for <f> and <p and the relation 
cf) = cf) hold only for this p = 1/4 case. Replacing exp(— 9) with exp# and interchanging p^ 
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with pn gives Pr. We get rid of P-(y) in these expressions by substituting P(y) — P+(y) 



for it; using (|2.16|) and the relation 



dy 1 lf„ cos 8k „ , . _ s 



27r cosh(y — 8) cosh y 4 J cosh 2 (^ 
giving us the remarkably simple relation 

2nP L {8) = ^le-° + [ dy P+{v) ; (2.20) 
2 7 cosh(y — 0) 

Pr(8) is the same with exp —8 replaced with exp6>. 

We can now find the free energy and hence the ground-state energy by minimizing 
the free energy subject to the constraints ( |2.16|) and ( |2.20|) [p!5| ^|,p| . This is standard by 
now; we provide a few more details in the appendix. The result is that the ground-state 
energy is 

E(R) = ~ j d8(e° ln(l + e'^) + e~ e ln(l + e" 6 ^)) (2.21) 
where the e a obey the integral equations 

fl(flH ^_|^_i_ ln(1 + e -M,, ) 

«< 9 > ^ "/ % t^V) (hi + «-"-') + m +e -">)) < 2 - 22 ' 

e«(6l) = e L (-(9). 

This set of integral equations is summarized by the diagram 

e e 1 e~ 9 

<s> — o — 

This is the result conjectured in The properties of these equations were analyzed in 
P,|I4]|, so our only comment is that the effective central charge c = 6RE(R)/n does indeed 



start at 4/5 in the UV and decreases smoothly to 7/10 in the IR. 
For general t, the conjecture of Jl[ is that the TBA system 

e e 1 2 t-4 e~ 

<s> — o — o o — <s> 

describes these flows. This does have all the required properties, so we expect that doing the 
TBA calculation with our S'-matrix will indeed give this result. Proving this is somewhat 
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tricky. Following G2], one can obtain equations for the eigenvalues and eigenvalue zeros. 
The problem is that the equations for the eigenvalue zeros have more solutions than are 
needed; some need to be thrown out. One can either do this by hand (as in |22| j) or by 
first studying the untruncated S'-matrix and then obtaining the RSOS result by adding 



appropriate chemical potentials to change the boundary conditions around the circle JT3|] • 
For the massive case, this enables the proof of the TBA equations for all t, but in the 
massless case, various technical complications arise with either method. It is possible to 
complete this proof, but it is complicated and we defer it to the future. However, in the 
infrared limit, Srl — ► 1 and the left and right particles decouple. The TBA calculation 



then is virtually identical to that of [Tj], and we obtain the correct IR central charge. We 
will do this calculation in the next section; we take it as a strong piece of evidence that 
our S'-matrices are correct for all t. 



3. The sine-Gordon model with imaginary potential 

In the previous paper 0, we discussed a number of properties of the sine-Gordon 
model 



S SG = J d 2 x [\{dcj>) 2 + y 2 cos (3 SG <i>\ 



at imaginary fi 2 . This describes the continuum limit of the low-temperature phase of the 
0(n) lattice model. It has the free-boson critical point at y 2 = 0, but it has another fixed 
point which in the 0(n) model corresponds to low temperature. In the phase between 
the two critical points, the excitations must be massless; these are the only ones which 
remain in the infrared limit. One way of studying this model is to analyze the ground-state 
energy on a cylinder E(R); in @] we discussed this by analytically continuing the high- 
temperature results and by lattice simulations. In an integrable model such as this, one 
calculates E{R) and the effective central charge c{MR) = 6RE(R)/n from the S'-matrix 
by the TBA. This model has central charges of 1 in both the ultraviolet limit and in the 
infrared limit |[23|| . 

The most plausible guess for the S'-matrix is given by the untruncated S'-matrix in 
and ( j2.4|) , where n and 0sg are related to t by 

n = 2 cos ( - J = 2 cos— — ; 3l r = , 3.1) 
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and where M oc | / u| ( - t+1 ^ 2 . This is for a number of reasons. The most compelling is that, 
as we showed in the previous section, this S-matrix can be reduced to give the minimal- 
model S-matrix; we expect this behavior from sine-Gordon at imaginary coupling, as was 
shown from large-t perturbative arguments in |24| (the exact symmetry to be used for this 
reduction may however not be the obvious one: this is a delicate issue which we discuss in 
sect. 4). 

Taking some simplifying limits gives further support for this S-matrix. When (j, — > 
(t — > oc and n = 2), the RL scattering is trivial. The L and R systems are decoupled for 
all M, not just in the IR limit. This system has an ordinary SU(2) symmetry, and was 
studied in ||. This system has c(MR) = 1 for all masses; this is what we expect from the 
1/t expansion done in 0. In the other extreme limit n = — 2 (// — > oo), the pathologies of 
the 0{— 2) model [0 lead one to expect that the S'-matrix description breaks down here, 
as indeed happens in our S'-matrix. As a final support for this S-matrix we now present 
computations at zero temperature but with a magnetic field, which turns out to be an 
easier problem. 

3.1. Adding a magnetic field 

We start with the sine-Gordon action coupled to an external U(l) gauge field 



S{A) 



J \{d^Y cos (3sg¥ d 2 x + J A^j^x. (3.2) 



Here j M is the U(l) current of the sine-Gordon model 



3p. = -z—tpvOutp 
normalized so that in the ordinary (/i 2 real) sine-Gordon model, 

n f a 0sg f°° dip 

Q = / Jodxi = — — / - — dxi (3.3) 
J 2tt dx x 

is the integer- valued soliton topological charge (the soliton (antisoliton) has then Q = 1 
(— 1)). Let us take A^ constant (i.e., space-time independent) and consider the correspond- 
ing specific vacuum energy £(A, fj 2 ) as a function of the "magnetic" field A = 

Before we turn to the scattering theory, it is worth looking at the action ( |3.2| ) from the 
perturbative (in /j, 2 ) point of view. Dimensional arguments as well as explicit perturbative 
calculations show that the magnetic field works as an infrared (IR) cutoff at scales ~ A 
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and therefore if A ^> fi 1+t the theory is in the ultraviolet regime. As a leading A — > oo 
approximation we set \i — in (1). Then a one-line calculation gives for the ground-state 
energy density 

£(A,0) = -^fA 2 . (3.4) 
In general it is convenient to define the effective central charge k(A, fi 2 ) as 

S(A,ii 2 ) = -— k(A,^ 2 ). (3.5) 

TV 

It is plain from the scaling considerations that k(A, fi 2 ) is actually a function of the di- 
mensionless variable 

£ = ^ 2 M 2/(1+t) 

and the perturbative series is in fact an expansion in £ 2 

oo 
1=0 

where in particular (t is related to the sine-Gordon parameter f3 as in ( |3.1| )) 

t 



ko 



k 2 = — 



1 + t 

1-t 



4 \i + tj r( — 

At this point let us suppose that the series (|3.6|) has some finite radius of convergence £o 
defining therefore an analytic function fc pe rt(£) at |^| < £o- In view of the subsequent Bethe 
ansatz analysis this suggestion seems plausible. We also call S peTt (A, fi 2 ) the ground state 
energy defined by this perturbative series 

A 2 

£ pe rt(A, ll 2 ) = fepert(0- ( 3 - 7 ) 

7T 

Note that this perturbative definition does not care if [i 2 is real or imaginary; it is inde- 
pendent of the physics of the corresponding field theories. 

We now discuss the scattering theory. Let us choose X\ as the space-like 2D coordinate 
and direct the constant field A^ along the "imaginary-time" coordinate Xq. The magnetic 
field perturbation in action (|3.2|) modifies the sine-Gordon Hamiltonian 7i s G to be 

H(A) = H sG - AQ (3.8) 
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Consider first the unitary massive sine-Gordon model {p 2 real in (|3.2|) ) and let m 
be the mass of the corresponding charged particle (soliton). As usual the on-mass-shell 
momenta (e,p) are parameterized in terms of rapidity 9 



e = m cosh 9 ; p = m sinh 9 

Under ( |3.8| ) every soliton (antisoliton) acquires additional energy A (—A). It is clear that 
if A > m the state without particles is no longer the ground state. The true vacuum 
contains a sea of positively-charged solitons which fill all possible states inside some (A- 
dependent) "Fermi interval" —B<9<B. There is non-trivial scattering between the 
solitons which certainly influences the structure of the ground-state sea. However, only 
one kind of particle participates and we can handle their factorized scattering as a diagonal 
one, the two-particle amplitude being —a(9) of (|2.5| ) with 

for the moment. Along the standard lines of the Bethe ansatz (BA) technique (see e.g. 
|2"5| - |2"7[] ) one arrives at the following expression for the specific ground state energy £^ Re \A) 
(the superscript Re is added to stress that currently we address the ordinary sine-Gordon 
with real coupling p, 2 ) 

/B 
cosh 9 p(9)d9. (3.9) 
-B 

The positive function p{9) (the particle density in the ground state) is defined inside the 
Fermi interval — B < 9 < B and solves the BA equation 

A - m cosh 9 = p(9)- [ <j>{9 - 9')p(9')d9' (3.10) 

J-B 

while the Fermi boundary B is determined by the boundary condition 

p(±B)=0. (3.11) 
The kernel <p(9) follows from the soliton-soliton scattering: 
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We are primarily interested in the UV (A — > oo) analytic structure of £( Re \A). Since 
B — > oo in this limit, it is convenient to transform (14) to the following Wiener- Hopf 



equation 28| , 25| - g7 



iAK + (0) + ime B K+(-i) 



to 



2 u — i 

for the meromorphic function v(u). Here 

K+(0) = y/2t/(l+t) 



e 2iuj ' B doj' 
uj + uj' K J K ! 2m 



(3.13) 




t \l + t 



t/2 



1 + t 



(3.14) 



and the kernel function a(u) reads explicitly 



,2iwA 



Y 



i(l + t)u 



r ( — ^ r + — ^ r ( i ^ 1 + t "> UJ 



(3.15) 



with 



A = -logt-— log(l + *). 



The integration contour C+ encircles all singularities on the positive imaginary axis. These 
are only the poles of a(u) at 

2in 



to 



1 + t ' 

The boundary condition ( |3.11|) now reads as 
iAK+(0) - l -^K + (-i) 



n = 0,1,2, 



(3.16) 



e a(w)w(w) — 



(3.17) 



while the ground state energy ( |3.9| ) becomes 



£ (i?e) (A) -£ (i?e) (0) 



me 



B 



2tv 



-K+(-i) x 



me"" . i / e ~"~- , ... da; 

AK + (0) - —K + (-i) + -a a,),,^) 

4 a; — z zm 



(3.18) 
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In the last equation the integration contour contains one extra pole at ui = i which one 
can pick up explicitly. We recognize its contribution \m 2 tan(7rt/2) as (with the opposite 
sign) the well-known bulk vacuum energy (without any magnetic field) of the sine-Gordon 
model [p7| , [29[1 , i.e. what we have denoted £^ Re \0) in eqn. 



£(*e) (0) = __ tan _ 



(3.19) 



The bulk term looks infinite for t an odd integer because in this case it depends on m 2 lnm 
|lj]. The remaining part of (|3.18|) , using the scaled definition ( |3.5|) , develops gular 
expansion in powers of (m/A) 4 ^ 1 ^ i.e. 



k( Re \A,m) =J2 K i (^M) 4V(1+t) 

1=0 

This series should be the same as the perturbative one ( |3.6|) provided that fl27j 



(3.20) 



l + t 



l + t 



l + t 



2/(1+*) 



(3.21) 



From a formal point of view one may interpret the BA system ( |3 . 9| ) — ( p7TT|) as a 
convenient way of summing up the perturbative expansion ( |3.6|) at real One would like a 
similar tool to sum up ( |3.6|) at £ purely imaginary, i.e. the alternating series. We show here 
that the massless scattering theory proposed above provides the necessary information to 
resum this series. Let us turn on the magnetic field perturbation (|3.8|) (say with A > 0) 
in our new world of massless charged particles. The energy spectrum 



e = p 



-P 



M 
2 

M 



for right-movers 
for left-movers 



is gapless and the positively charged particles (say the u ones) are always excited in the 
ground state once A > 0. Now the right- and left-movers fill respectively the semi-infinite 
Fermi intervals — oo < 9 < B and — B < 6 < oo with some Fermi boundary B ~ logA/M. 
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Again it is a straightforward Bethe ansatz exercise to derive the following system of integral 
equations 

B oo 

Me 6 f f 

qa- — = PR {d) -J mo- o') PR {d')de' - J mo - o') PL {e')de'- 

(3.22) 

oo B v ' 

Me~ 6 f f 

QA — — = PL (e)- J Mo-o') PL {e')de' - J mo -o') PR {o')d9'. 

-B -oo 

The positive functions P r(0) and P l{0) are defined in the Fermi intervals oo < 9 < B and 
—B < 9 < oo respectively, and are restricted by the boundary conditions 

PR (B) =p L (-B) =0 (3.23) 



In eqn. (|3.22|) we introduced the U(l) charges ±Q of the massless u and d particles, which 
we see no reason to fix in advance. The kernels <pi(9) and <p2 {9) are related to the LL and 
RL diagonal uu — > uu scattering amplitudes respectively (see ( |2.5| ) and ( |2.4| ) where now 



M = !/(*-!)) 



Id, f e- e sinh^^ du 



2cosh^sinh^i^ 2tt 



^0) = ^—^loga RL (9) = - [ 
2m d9 J 



Id, _ r e^sinh^f du 

2 cosh if sinh n{t ~ 1)uJ 2tt ' 



(3.24) 



Finally, the ground state energy, which we now call £^ Im \A), is evaluated as follows 



M r B 

£ {Im) (A) - £V m \0) = -— e 9 PR (9)d9 7 (3.25) 

27T J-oo 

where we have taken into account the obvious symmetry Pr (9) = P l(—9). 

Although this new BA system (|3.22|) -( |3.25| ) has a rather different form from that of 
eqns. ( |3.9| )-( |3T2"0|) , it is easy to relate the two in the UV region A — * oo where in both 
systems B — > oo. In this limit the right and left Fermi intervals have a broad overlap at 
—B < 9 < B. Near say the right Fermi boundary 9 ~ B (where the main contribution 
to ( |3.25| ) comes from) we can forget about the left one and solve ( [3.22 ) for P l(9) by the 



Fourier transform with B — > oo. The resulting equation for Pr {9) is 

Me f B 
rQA-—- = PR {9) - $(9- 9') PR {9')d9' (3.26) 

^ J-oo 
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where in terms of the Fourier transforms 

[&(*>)] a 



<&(w) = (pi(uj) + 



sinh 



7r(t-l)cj 



l-0i(w) 2 cosh ^ sinh ^ 



(3.27) 



(compare with eqn. fl3.12f) ) and 



r = 1 + 



«(0) 



t - 1 



l-0l(O) t 

It coincides precisely with the corresponding limit B — > oo of eq.(14) provided 



t- 1 

M = m 



(3.28) 



For more detailed information we again transform the system ( |3.22|) - (|3.25|) into 
Wiener-Hopf form: 



u{u) 



iAK + (0) iMe B K+(-i) 
-rQ ±^ + — + 



to 



u — I 



duo' 



-a(u')u(u>')- 

c uj + uj' v J K J 2m 



(3.29) 



with the same K+(0) and K+(—i) as in eqn. (|3.14j) and the same integration contour C+. 
Likewise the boundary condition ( |3.23| ) turns to 



iMe B 

irQAK+{0) ^—K + (-i) 



e 2lU}B a(oo)u(oo 



duo 
' 2tt? 



(3.30) 



while the ground state energy reads 



£( Im \A)-£( Im \0) 



Me 



B 



2tt 

Me B 



K+(-i) x 



Me B f 
rQAK + (0) - —^—K + (-i) + / 



e 2iu,B Jw 

a(u;)'u(a0 

o , u - i 2m 



(3.31) 



With the relations ( |3.28j ) the only difference between the two systems ( |3.13| )-( |3.18[ ) and 
( p.29| )-( [3.31D is in the kernel function which is now a bit modified 



a(uj) 



sinh^ , 

7— a(uj) 

sinh ^ V ; 



(3.32) 



where a(u) is that of eqn. ( |3.15| ). Picking up the pole u = i in ( |3.31| ) we find the bulk 
vacuum energy of the sine-Gordon model at imaginary /j, 2 : 

M 2 



4cos(7rt/2) 
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At t integer this is the same bulk term as in the corresponding (reduced) interpolating 
flows between the minimal models QTJ] (as speculated in p5|), a result we used in our first 
paper. 



The analytic structure of the remaining part of eqn. ( |3.31| ), i.e. of 



7T 



is quite regular and controlled by the poles of the kernel function a{ui). These form just 
the same set ( |3.16| ) as that of a{u) (note that the poles of the extra multiplier in ( |3.32| ) 
are all cancelled by the appropriate zeros in ( |3.15[ )). Moreover, the residues are almost the 
same, only the sign alternating at n odd in (|3.17| ). This obviously results in the following 
analytic structure 

00 / M \ 4Z/(l+t) 

k^(AM) = J2(-) l K l ( " 

1— n \ 



-Y 



x rQA J 
1=0 \ ^ / 

with precisely the same coefficients K\ as in expansion ( |3.20|) . 

We conclude that up to the known bulk vacuum energy contributions the massless 
BA system ( |3.22j )-( p.25| ) gives the correct analytic continuation of the massive one (|3~^|)- 
( P-12| ) to purely imaginary £, providing (|3.28| ) holds. In particular, the low-temperature 
mass scale M is in the same relation ( |3.21| ) with the absolute value of ^i 2 . 

The success of this computation justifies the t — > t — 1 shift, the LR doubling and the 
shift of 9 in the LR matrix elements. It seems to make the above S'-matrix the right choice 
for the non-unitary sine-Gordon model. However, as we shall now see, the TBA does not 
lead to the correct form of the ground state energy on a cylinder (i.e. the finite-temperature 
free energy). 

3.2. The case t = 2 

The S'-matrix appears natural in the study of the n=0 point too. At the UV fixed 
point, this model has N=2 supersymmetry, but in the low-temperature phase, the su- 
persymmetry is spontaneously broken (this is possible thanks to the non-unitarity) . The 
spectrum should consist of four massless Goldstinos, one for each supersymmetry. The 
left movers and the right movers each form a (it, d) doublet. At the zero-temperature 
fixed-point, this model has c = —2, so in the infrared limit, the Goldstinos become the 
well-known ?y£ ghost system. This is a non-unitary version of what happens in the flow 
between the tricritical Ising model and the Ising model []3(| (the t = 4 case in the previous 
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section), where a Goldstino resulting from the spontaneously-broken N=l supersymmetry 
becomes the Ising-model free-fermion in the infrared limit. One can write down a Landau- 
Ginzburg action for the N=2 case. The N=2 superfield has two fermions and two bosons; 
to obtain the Goldstino interactions one can integrate out the bosons, as was done in the 



N=l case in ||30| . One also expects the LL and RR S'-matrices to be the identity, because 
the LL scattering essentially describes only IR properties, and as we mentioned above, the 
n = system in this limit reduces to the free ?y£ system. Srl diagonal then follows from 
factorizability. Our proposed S'-matrix has these properties for fi = 1. 

We can now consider the scattering theory at finite temperature and, via the TBA, 
extract the ground state energy on a cylinder. Surprisingly, it does not have the expected 
behavior. We provide the TBA equations for this diagonal scattering in the appendix. 
Although c e g = 1 in the IR, it is equal to .6 in the UV. 

Using this S-matrix, we have also calculated tr(— 1) F in the t = 2 case. This corre- 
sponds to changing the boundary conditions around the cylinder from the Neveu-Schwarz 
sector to Ramond. In a supersymmetric theory, lntr(— 1) F does not depend on L, so it 
should vanish in a TBA calculation. We do this calculation by inserting chemical potentials 
of m/R for the Goldstinos, and in the infrared limit it does not vanish but instead gives 
c = — 2. This is not the ground-state energy, but in fact is the actual central charge of the 
IR fixed point. (The effective c = 1 results from subtracting 1/12 times the dimension of 
the lowest-dimension operator in this non- unitary theory.) When we move away from the 
IR limit, we find striking behavior: this quantity has a discontinuity at mR ~ 3.6 (near to 
r2 of the previous paper), where it takes value c ~ —4.10. For mR smaller the TBA does 
not seem to make any sense except when postulating that there are no particles left, which 
corresponds to an effective central charge equal to zero. This strange behavior is actually 
identical to what was observed in lattice simulations, where the singularity occurs due to 
the level crossing. Compare in particular with figure 5 of the preceding paper. Thus it 
seems that the S'-matrix leads to the correct ground state energy in the R sector, but not 
in the NS sector. 



3.3. General t 

To find the ground-state energy using the TBA for general t, we must diagonalize the 
"transfer matrix" as before. This time, we use the algebraic Bethe ansatz ||31|1 - Just as 
we did in the minimal-model case, we can reduce this to an all-LL problem by shifting 
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the rapidities of all the right-movers by iTc/2/i. This is the case studied in [31], and the 
eigenvalues Al(8, {Oi}) and A R (9, {9i}) are 



(3.33) 



Nl Nr m , _ n\ 

Al(#; {eh M) =Ua(e- of) *rl{o - *?) ]J ^— 1| 

3=1 k=l r=l ^ r ' 

Nl Nr m ( 0\ 

A R (9; {6}, {y}) = j] a RL (9 - Of) ]J a(9 - 9*)n f[ aRL f r _ " 

j=l k=l r=l RL ^ r > 

where m is an integer < m < Nl + Nr, and the {y r } are solutions of 

f| b (Vr ~ 9 j) t4 b RL (y r -6£) = yr a(yr ~ ViMm ~ Vr) (3 M) 

f \ a(y r - 9}) 11 a RL (y r - 6%) f = \ a(y t - y r )b(y r - Vi ) 1 ' ' 

One sees that these equations are invariant under exchange of left and right by making the 
substitution y r — > y r + in/ (2 fx). 

As in the last section, one defines particle densities Pl(9) and p R (9) for the left- and 
right-movers. For general fi, there are many more types of solutions of ( |3.34|) than the 



simple ones for its analog ( 2.15 ); these are called string solutions [32]. One then defines 
a density for each type of string solution and then proceeds as before. There is, however, 
a major complication in our case which makes the TBA somewhat more intricate. We 
discuss this in the appendix. We were only able to find a result for the t = 3 case, and 
while it has c e g = 1 in the IR, it is .2 in the UV. Thus this seems to be unsuccessful like 
the t = 2 case. 

We can, however, analyze these equations successfully in the IR (M — > oo) limit for 
any integer > 2, and derive that ci R = 1 for the untruncated case, and cj R = 1 — 6/t(t — 1) 
for the RSOS models. In this limit, S R l = 1 and the left and right sectors decouple. The 
TBA for each of these sectors is essentially the same as that for the massive sine-Gordon 



models discussed in [0, and is described by the diagram 

e -e L 

<g>- -o- -o- 



for the left, with the same for the right. There are open dots in the diagram. This 
system has ci R = 1. 
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This result allows us to derive cm for the RSOS S-matrices discussed in the previous 
section. These RSOS S'-matrices were obtained from the untruncated one, and the TBA 
calculations are identical except that appropriate twisted boundary conditions must be 
placed around the cylinder to reproduce the RSOS model (this is the analog of the charge 



at infinity in conformal models). As discussed in [fT3| , this has the effect of removing the 
three nodes at the right end of the diagram. Thus the above diagram becomes the one at 
the end of the last section, but with the right-moving particle removed. This is of course 
the conjectured answer for the flows between minimal models, and is a good consistency 
check on the S'-matrices discussed in the last section. 



4. Discussion 

One thing we see is clear. At any coupling in sine- Gordon we have c = 1 in the IR 
limit, but at the values t = 2, 3 where we can do the TBA calculation, c 7^ 1 in the UV 
(for the trivial case t — > 00 we still get cjjv — !)• Thus, at least for t = 2, 3 the S-matrix 
cannot be used to describe, via the TBA, the flow of the running central charge between 
the sine-Gordon fixed points, which have c = 1 in the IR and UV. 

A first possibility is that the S-matrix is incorrect for every value of t. However 
we know that it gives the right TBA after truncation. This seems compatible only if 
additional massive particles are present, but disappear in the truncation. For this to 
happen, the particles should have a non-vanishing charge, and thus affect the magnetic 
computation. Since we have found this perfectly consistent, the existence of additional 
massive particles is doubtful. Similarly, the computations in the magnetic and truncated 
cases seem to exclude the presence of additional CDD factors. 

A second possibility is that the S-matrix is incorrect only for t < 3, which is unfortu- 
nately the only range we were able to study with the TBA. One hint towards this possibility 
is that for t < 3, there is a pole in the S-matrix which in the massive case would correspond 
to another state in the spectrum. For example, for t = 2, where Srl = ±tanh(#/2 + z7r/4), 
this pole is at 9 = in/2. In this massless case it is conceivable that it corresponds to the 
massive sine-Gordon particle, which now is an unstable resonance. Thus there may be 
some sort of phase transition at t = 3, where this pole crosses the "physical" threshold of 
9 = in. 

A variant would be that the S-matrix is correct for t < 3 but cannot be used to get the 
running central charge via the TBA. This is conceivable because the S-matrix is defined 
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only in an infinite space-time volume whereas the TBA applies at finite Euclidan time 
(except in the IR limit, where it gives the expected answer). Let us discuss this further. 
Consider first t = 2. Recall that N=2 supersymmetry is not explicitly broken in the flow. 
One therefore could expect that it is realized nonlinear ly in the scattering theory. However, 
the fermion number is unambiguously defined along the entire flow, including at the IR 
fixed point, and is proportional to to the topological charge (this can be checked explicitly 
for various L-leg polymer observables in the dilute and dense cases). For sine-Gordon 
solitons in the massive flow one has F = ±1/2. For rj,£ and their complex conjugates 
in the dense N=2 case one has F = ±1. The value of F for the massless particles in 
the S-matrix is therefore determined, and F should act linearly on multiparticle states. 
Therefore one has unambiguously that the Witten index tr(— 1) F (which gives Ramond 
boundary conditions) is expressed in the scattering theory for t = 2 by adding fugacities 
A = — 1 for the fermions. But it is easy to check that with such boundary conditions in 
space and time directions, there is no cancellation of various contributions in the scattering 
theory, and for any value of MR larger than the corresponding partition function is not 
the expected constant, as discussed at the end of sect. 3.2. This suggests that N = 2 
supersymmetry cannot be realized at all in our scattering theory. Such a possibility may 
not be surprising. Recall that even though levels cross the ground state of the Ramond 
sector and spontaneously break the N = 2 supersymmetry, the Witten index remains 
equal to 2 along the flow. We set it by hand to zero at the IR fixed point since the bulk 
term is there infinite and the levels contributing to this index are infinitely far from the 
ground state. Reaching this IR fixed point is thus a rather singular process, and since 
the scattering theory is a perturbation of this IR fixed point, it appears reasonable that it 
cannot be used "to go all the way" to the ultraviolet model. More precisely we gave in the 
previous paper numerical indications that in the Neveu-Schwarz sector the ground-state 
energy has singularities at r2 and r±, and maybe at all values of r2 n as well. If this is true, 
one may be in a phase that cannot be described by the TBA and scattering theory. Such 
a situation is expected when the virial expansion does not converge |33|] . 

The case t = 3 is a little marginal. The associated minimal model is in fact massive 
- it is the Ising model flowing from c = 1/2 to the trivial c = fixed point. One does does 
not expect to find an S'-matrix describing a perturbation of the IR fixed point to be able 
to describe the flow, because it has only trivial degrees of freedom on which to act. This 
is reproduced by the reduction of our sine-Gordon S'-matrix, which gives a model with 
no particles and thus a ground state energy always equal to zero (and not reproducing 
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cjjv = 1/2). As an aside notice that S'-matrices with U q isl{2) quantum group symmetries 
acting on spin 1/2 representations can be interpreted directly, using another representation 
of the Temper ley-Lieb algebra, as ^-matrices describing scattering in an 0(n') symmetric 
theory with n' analytically continued to n' = —q' — q'~ x (for the massive case, see |34|,|35|] ) . 



For t = 3 we have q' = z, which gives n' = 0. The result fits with a naive Goldstone picture, 
where in the low-temperature phase of the 0{n = 1) model only degrees of freedom with 
0(n' = n — 1 = 0) symmetry are present. 

At the present time we can only hope that our S'-matrix gives the correct finite- 
temperature results for t > 3. 

To conclude we discuss briefly the strange symmetry pattern of the problem, which 
may be the origin of this confusion. We expect the minimal model to be obtained by 
twisting and reduction of the full sine-Gordon theory, using the U q sl(2) quantum symmetry 
with q = — exp(— in/t) (the subalgebra U q sl{2) with charges of zero dimension in the 
twisted theory appears on the lattice) . For the scattering theory, this symmetry is expected 
to act nonlinearly. On the other hand, the minimal S'-matrices we have conjectured above 
are obviously related to U q >sl(2) symmetry, where q' = — exp(— in/(t — 1)), acting linearly. 
"Undoing" the truncation as we did is equivalent to assuming that the linear U q >sl(2) 
symmetry is a nonlinear realization of the subalgebra U q sl{2). In that case the scattering 
theory for the full sine-Gordon is indeed obtained by going from RSOS to vertex basis. 
This certainly makes sense from the point of view of the topological charge, i.e. the Cartan 
subalgebra. Indeed, we showed above that at the IR fixed point, the central charge of 
the minimal model was obtained from the sine-Gordon TBA by introducing fugacities 
A = — q f±1 . On the other hand, in the lattice model as well as for sine-Gordon in the UV, 
the central charge of minimal models is obtained by taking traces of evolution operators 
with the insertion of <p where Q is the topological charge ( |3.4| ) 



2n 



dx d 7 



This is of course the soliton number in sine-Gordon. It is as well proportional to the 
magnetic number in the Coulomb gas and solid-on-solid models, the spin in the 0(n) 
vertex model, and the Cartan generator H(S Z ) of U q sl{2). From this one expects 
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for the topological charge carried by the massless particles of our scattering theory. We 
found this result independently in the magnetic calculation, so everything seems consistent. 
The nonlinear realization of the other generators remains however to be found. It may be 
that the mechanism of truncation in this massless case proceeds very differently from our 
naive approach. In such a case, the way we have "undone" the quantum-group truncation 
may be incorrect, and the non-unitary sine-Gordon S'-matrix very different. 
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Appendix A. The TBA equations 

Starting with equations for the level densities of the form 



2nP a {6) = lu a {9) +Y.f dd> 'Pb(0')<t>ab(0 - 9'), (A.l) 



b 

one then minimizes the free energy, and obtains the ground-state energy 



E X (R) = -J2^ J ddv a (d)ln(l + \ a e-^), (A.2) 



where the functions e a obey the TBA equations [15 



^ b (e-e')Hi + x b e-^^), (A.3) 

and where the A a are fugacities. It is important to notice that the kernel which appears 
in ( |A.3|) is the transpose of that which appears in ( |A.1|) ; i.e. 4>Jj = 4>ji. In most cases, the 
kernel is symmetric, but we will see below that in our sine- Gordon case it is not. To find 
the central charge in the UV limit, one can rewrite ( |A.2|) in terms of a sum of dilogarithms 
| |15| |; this argument only works when <fi = <p T . 

In the simple example /i = 1 from section 2, the scattering is diagonal so (p a b = 
— id In S a b(9)/d9. This leaves us <Prl = —1/ cosh 6*, and <Pll = 4>RR = 0. There are 
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two right-moving particles with vr(6) = (M/2)exp#, and two left-movers with = 
(M/2) exp(— 0). The analysis of dilogarithms then gives cjjy = .6 and cm = 1. To 
calculate tr(— 1) F , one includes fugacities of —1 for all the particles. In the IR limit, this 
gives c = — 2; the UV limit is ambiguous because of the discontinuity of the solution at 
mR ~ 3.6. 

We would now like to study the TBA equations for general p. There is a unusual 
complication here. Let us study this in the simplest example p = 1/2, where the left- 
hand-side of (pPID is just (-1)^+^* = l. @ Then for large N L + N R , the solutions of 
( |3.34| ) are of the form y = z + z7r/2, or y = Zq — iiv/2, where zq and Zq are real. The 
derivative of the log of ( |3.34| ) then yields 

±2^p t {z % ) = I de p " {e) - - [ de Pf ] (A.4) 



cosh(zi — 6) J cosh(^j — 6) 

where i = or 0. The ± in front is normally fixed by demanding that Pj be positive. 
This is where the complication arises: the sign of the right-hand-side depends on the value 
Zi. This is easy to see, using the fact that by symmetry pl(0) = pr(—8). Thus we must 
replace the ± by sgn(zi). The TBA relations can then be derived, but many features are 
more complicated. The kernels a 6 involve this irritating sgn, and are not symmetric in 
o and b. As a result, the densities one derives have discontinuities. However, the TBA 
analysis can still be formally done, although we are not convinced that it means anything. 
In particular, when one tries to add fugacties to truncate to the minimal models, the 
equations produce a complex answer even though the energy should remain real. In any 
case, we have studied the TBA equations for this case p = 1/2. We obtain an effective 
central charge which starts at .2 in the UV and increases smoothly and monotonically to 1 
in the IR. The tricks one usually uses to rewrite the ground-state energy in the ultraviolet 
limit as a sum of dilogarithms must be modified because of the discontinuities, but they 
give the result of .2. 

In the general case p ^ \ outside of the IR limit, things are even more complicated. 
The problems with signs are worse, and it is not even clear how to classify the solutions 
of the equation ( [2.15| ). The usual string hypothesis [^2J is not valid, because it destroys 
the left-right symmetry. We have made a modified string hypothesis for 1/p an integer, 
where there is an antistring for every string. This restores the left-right symmetry. We 
have studied the resulting TBA systems, but not found anything of interest yet. Again, 



See [36|, sec 6.1, for a similar calculation done in detail. 
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it is not clear whether this has any meaning. The fact that these extra antistrings must 
disappear in the IR limit where we understand the system in fact leads us to believe that 
the solutions of these TBA equations may not be simply described in terms of strings. 
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